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INTRODUCTION 


Over the past three years we have been developing a finite element based procedure for 
the solution of high speed viscous compressible flows. The approach followed has been to 
compute steady state solutions via a false transient, using an explicit time stepping 
scheme, and to attempt to improve the solution quality by incorporating adaptive mesh 
procedures. This strategy has been shown to be very successful in the context of the 
compressible Euler equations and also for some simple compressible viscous flows. The 
main thrust of our work over the past twelve months has been to work on the extension of 
the approach to the solution of some realistic compressible viscous flows. When flows at 
high Reynolds number are investigated, it soon becomes apparent that explicit techniques 
have to be supplemented if they are to deal effectively with the large variations in element 
size and aspect ratio which characterise the computational grids necessary for adequate 
resolution of the primary flow features. For this reason, our basic Taylor -Galerkin 
solution algorithm has been re-written in an expicit/implicit form. A structured grid is 
laid down in the vicinity of solid walls and the remainder of the flow field is covered by 
an unstructured assembly of triangles. In the unstructured grid region the basic explicit 
solver is employed, while the implicit version is used in the structured portion of the 
grid. The solution of the implicit equation system is achieved by directional splitting, 
using a block tridiagonal solver. Using this approach, solutions have been computed for 
the problems of (i) flow past a flat plate, M-3, Re- 1000 (ii) shock/ boundary layer 
interaction, M-2, Re-296000 (ill) flow over a compression corner, M- 11.68, Re-246000 
Uv) uniform flow past a circular cylinder, M-6.34, Re-39770. A summary of the results 
produced is included in this report and demonstrates the numerical performance of the 
scheme. 


IMPLICIT/EXPLICIT SOLUTION SCHEME 


A description of the scheme will be given for the Euler equations. The extension [1,2] to 
the Navier-Stokes equations follows directly, but the algebra becomes rather messy. The 
form of the scheme was suggested by the original work of Lerat et al [3,4]. The governing 
equations are assumed in the form 
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where, in the usual notation, U denotes the vector of conserved variables and Fj denotes 
the flux component in direction ij of a Cartesian coordinate system. Using a Taylor series 
expansion about time t - t„ gives 
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and, using the governing equation, we can write 




vg 

= 






The time stepping scheme then becomes 
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The solution of this equation system is achieved in a two-stage manner: 


STAGE 1... explicit predictor 
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This stage is performed using the two step Taylor-Galerkin code that has formed the basis 
of our previous work. 

STAGE 2. .implicit corrector 

The implicit corrector is only employed on the structured grid and the equation system is 
solved by directional splitting in the form 
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To reduce the computational cost, these equations are formed using the finite difference 
approach. Note however that at steady-state the solution satisfies 

= O 

APPLICATIONS 

The results produced by applying the above algorithm to a series of viscous flow problems 
will now be shown. 

(i) flow over a flat plate, M-3, Re- 1000. A regular mesh of 56*50 points and 5600 elements 
was employed ( see Figure 1 ). Points were clustered near the leading edge and near the 
plate. Isothermal wall conditions were applied on the plate with symmetry conditions 
ahead of the nose of the plate. The use of a regular mesh means that the implicit scheme 
can be applied everywhere or only locally. Figures 2 and 3 show the density and pressure 
contours when 15 layers near the plate are treated implicitly with a Courant number C-5. 
A comparison between the behaviour of the fully implicit, mixed implicit/explicit and 
fully explicit schemes is made in Figures 4-7 in terms of the convergence behavior, eiit 
density and velocity profiles and plate pressure distribution. Figures 8-10 show the 
comparison between the results of Carter 151 and the computations performed at ALB, LaRC 
for the same quantities. It may be observed that the rates of convergence of the fully 
implicit and the mixed explicit/implicit schemes are approximately equal, while the fully 
explicit algorithm shows a much slower convergence behaviour with peculiar oscillatory 
behaviour ( which may be caused by an inappropriate artificial viscosity model for highly 
stretched grids ). 

(ii) shock/boundary layer interaction, M-2, Re-296000. The computation corresponds to 
the experiments of Hakkinen et al (61 and the Reynolds number is based upon the distance 
of the shock impingement point from the leading edge of the plate. The shock is of 
sufficient strength to cause separation of the developing laminar boundary layer. The first 
mesh employed consists of a regular grid of 32*45 points, equally spaced in the x^ 
direction and exponentially stretched in the x^ direction. The grid is shown in Figure 1 1 
and the computed pressure and Mach number contours produced by the fully implicit 
procedure are shown in Figures 12 and 13- Adaptive remeshing was now performed, 
producing the mesh of 3293 points shown in Figure 14. On this grid the solution was 
advanced at C-5 in the 21 layers nearest the plate and the final pressure and Mach number 
contours are given in Figures 15 and 16. The computed variation in the static pressure and 
the skin friction along the plate are shown in Figures 17 and 18 respectively. The results 
obtained on the second mesh have been found to agree well with the computed results of 
MacCormack [71 and Dawes |8J. Details of the velocity profiles in the vicinity of the 
separation point can be seen in Figure 19. 

(iii) flow over a compression corner, M-11.68, Re-246000. The intention was to compute 
the configuration tested by Holden (9J but it has been recently discovered that the 
computed compression corner angle was 1 4.5 degrees which differed from the 1 5 degrees 
in the actual experiment. Consequently, this problem is curently being recomputed for the 



correct geometry and with an Improved mesh. We present the results for the incorrect 
geometry as an eiample of the computational performance. An initial solution was 
produced on a regular mesh, equally spaced in the x ^ direction and exponentiaily 
stretched in the x x direction. The computed solution on this mesh was used to start the 
adaptive refinement process. The final mesh employed is shown in Figure 20 and the 
solution was advanced with C-5 in the 24 structured layers along the wail. Density and 
pressure contours are shown in Figures 21 and 22 and the distributions of pressure and 
skin friction along the wall are given in Figures 23 and 24. The computed region of 
separated flow is smaller than that observed experimentally which suggests insufficient 
mesh resolution in the vicinity of the corner. 

(iv) undisturbed flow past a cylinder, M-6.34, Re-39770. The problem was initially solved 
inviscidly and the mesh adapted for the bow shock. A structured mesh was then 
constructed In the vicinity of the cylinder and consisted of 15 layers of elements. The 
final mesh employed had 4312 points. The solution was advanced with C-5 in the 
structured grid region and explicitly over the remainder of the grid. The final mesh and 
the computed contours of density and Mach number are shown in Figures 25-27. The 
computed distributions of pressure and heating rate over the cylinder surface are 
displayed in Figures 28 and 29 and are in good agreement with the results produced for 
this problem at ALB, LaRC. 

CONCLUSIONS 

A mixed explicit/implicit algorithm for the solution of problems involving high speed 
compressible viscous flows has been presented. Initial tests with the method have proved 
to be encouraging and the approach is being tested and extended further. In particular, 
the need for employing a structured grid and the finite difference discretisation in the 
implicit region is being reviewed with the objective of producing an implicit method for 
unstructured grids. 
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